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ABSTRACT 

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by 
Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of central 
galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass 
assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star 
formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula 
that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass 
and redshift. We use this model to make predictions for various aspects of the star formation rates 
of central galaxies across cosmic time. Our main findings are the following. (1) The specific star 
formation rate (SSFR) at high z increases rapidly with increasing redshift [oc (1 + z)^-^] for halos of 
a given mass and only slowly with halo mass (cx M^-^"^) at a given z, in almost perfect agreement 
with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation 
rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at 
a constant value, ~ 10~^'^/).^yr~^, independent of halo mass or the final stellar mass of the galaxy. 
However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of 
the stars in the present-day Universe were formed in halos with lO"'^^'^ /i~^Mq < Mh < lO"'^^''^ /i~^M0 
in the redshift range 0.4 < z < 1.9. (4) The star formation efficiencies (SEE) of central galaxies 
reveal a 'downsizing' behavior, in that the halo 'quenching' mass, at which the SFE peaks, shifts from 
~ 1012-5 /j-^Mq at z > 3.5 to - lO"-^ h~^MQ at z = 0. (5) At redshift z > 2.5 more than 99% of 
the stars in the progenitors of massive galaxies are formed in situ, and this fraction decreases as a 
function of redshift, becoming ^ 60% at z = 0. For a Milky Way sized halo of ~ 10^^ h~^MQ 
more than 80% of all the stars are formed in situ, as opposed to having been accreted from satellite 
galaxies. 

Subject headings: cosmology: dark matter - galaxies: formation - galaxies: halos 



1. INTRODUCTION 

Recent years have seen dramatic progress in estab- 
lishing the connection between galaxies and dark mat- 
ter halos, as parameterized either via the conditional 
luminosity function (CLE) (e.g., Yang et al. 2003) or 
the halo occupation distribution (HOD) (e.g., Jing et 
al. 1998; Peacock & Smith 2000). This galaxy-dark 
matter connection describes how galaxies with differ- 
ent properties occupy halos of different mass, and con- 
tains important information about how galaxies form and 
evolve in dark matter halos. In practice, the methods 
that have been used to constrain the galaxy-dark mat- 
ter connection (galaxy clustering, galaxy-galaxy lensing, 
galaxy group catalogues, abundance matching, satellite 
kinematics) use the fact that the halo properties, such 
as mass function, mass profile, and clustering, are well 
understood in the current ACDM model of structure for- 
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mation (see Mo, et al. 2010 for a concise review). 

At low redshift, large redshift surveys, such as the two- 
degree Field Galaxy Redshift Survey (2dFGRS; CoUess 
et al. 2001) and the Sloan Digital Sky Survey (York et al. 
2000) have been used to establish reliable links regard- 
ing how galaxies with different properties are distributed 
in halos of different masses (e.g. Jing et al. 1998; Pea- 
cock & Smith 2000; Berlind & Weinberg 2002; Yang et 
al. 2003; van den Bosch et al. 2003, 2007; Zheng et 
al. 2005; Tinker et al. 2005; Mandelbaum et al. 2006; 
Brown et al. 2008; More et al. 2009, 2011; Cacciato 
ct al. 2009, 2013; Neistein et al. 2011a,b; Avila-Reese 
& Firmani 2011). At intermediate redshift, z ^ 1, rela- 
tively large redshift surveys, such as the DEEP2 survey 
(Davis et al. 2003), the COMBO-17 survey (Wolf et al. 
2004), VVDS (Le Fevre et al. 2005), and zCOSMOS 
(Lilly et al. 2007), have also prompted a series of inves- 
tigations into the galaxy-dark matter connection and its 
evolution between z ~ 1 and the present (e.g.. Bullock et 
al. 2002; Moustakas & Somerville 2002; Yan et al. 2003; 
Zheng 2004; Lee et al. 2006; Hamana et al. 2006; Cooray 
2005, 2006; Cooray & Ouchi 2006; Conroy et al. 2005, 
2007; White et al. 2007; Zheng et al. 2007; Conroy & 
Wechslcr 2009; Wang & Jing 2010; Wetzel & White 2010; 
Wake et al. 2011; Leauthaud et al. 2012). However, at 
higher redshifts, especially beyond z ~ 2, reliable clus- 
tering measurements are not available, and the data is 
limited to estimates of the luminosity / stellar mass func- 
tions of galaxies (e.g. Drory et al. 2005; Fontana et 
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al. 2006; Perez- Gonzalez et al. 2008; Marchesini et al. 
2009; Stark et al. 2009; Bouwens et al. 2011), often 
with large discrepancies among different measurements 
(see Marchesini et al. 2009). It is thus not possible to 
carry out the same kind of HOD / CLF analyses for high-z 
galaxies as for galaxies at low z. Nevertheless, attempts 
have been made to establish the relation between galax- 
ies and their dark matter halos out to high z using a 
technique known as abundance matching (introduced by 
Mo & Fukugita 1996 and Mo, Mao & White 1999), in 
which galaxies of a given luminosity (or stellar mass) are 
linked to dark matter halos of a given mass by match- 
ing the observed abundance of the galaxies to the halo 
abundance obtained from the halo mass function (typi- 
cally also accounting for subhalos) (e.g., Vale & Ostrikcr 
2004, 2006; Conroy et al. 2006; Shankar et al. 2006; 
Conroy & Wechslcr 2009; Moster et al. 2010; Guo et al. 
2010; Behroozi et al. 2010). 

However, as pointed out in Yang et al. (2012; here- 
after Y12), although subhalo abundance matching yields 
galaxy correlation functions that are in remarkably good 
agreement with observations (e.g., Conroy et al. 2006; 
Guo et al. 2010; Wang & Jing 2010), it suffers from 
the following two problems. First, assuming that the 
stellar masses of satellite galaxies depend only on their 
halo mass at accretion implies either that the relation 
between central galaxy and halo mass is independent of 
the time when the subhalo is accreted, or that the ef- 
fects of different accretion times and subsequent evolu- 
tions in different hosts conspire to give a stellar mass 
that depends only on the mass of the subhalo at accre- 
tion. Second, as the satellite galaxies are forced to be 
linked with the subhalos that survive in simulations, no 
satellite galaxies are allowed to be associated with subha- 
los that have been disrupted in the N-body simulations. 
To circumvent these inconsistencies, Y12 proposed a self- 
consistent model that properly included the fact that (1) 
the relation between stellar mass and halo mass of cen- 
tral galaxies depends on z\ (2) the properties of satellite 
galaxies depend not only on the host halo mass at accre- 
tion, ma, but also on the accretion redshift, z^; and (3) 
after accretion a satellite galaxy may lose or gain stellar 
mass and even be disrupted due to tidal stripping and 
disruption. Based on the host halo and subhalo accretion 
models provided in Zhao et al. (2009; see also Zhao et 
al. 2003) and Yang et al. (2011), Y12 obtained the con- 
ditional stellar mass functions (CSMFs) for both central 
and satellite galaxies as functions of redshift assuming 
a number of popular ACDM cosmologies. In particular, 
the mass assembly histories of central galaxies, the pop- 
ulation of accreted satellite galaxies, and the fraction of 
surviving satellite galaxies are all well constrained. 

With the results obtained in Y12, it is straightforward 
to obtain the star formation histories (hereafter SFHs) 
of galaxies, especially for the central galaxies, in halos of 
different masses at different redshifts. Indeed, along this 
line, a couple of recent investigations have tried to model 
the SFHs of galaxies in different halos using a subhalo 
abundance matching (SHAM) method (e.g., Moster et 
al. 2013; Behroozi et al. 2012, 2013; Wang et al. 2012). 
In all these investigations, assumptions have to be made 
about the satellite evolution, their contribution to the 
masses of central galaxies, as well as the intra-halo stars. 
In this paper, we present our own modeling of the SFHs 



of (central) galaxies using a self-consistent model that is 
not hampered by the shortcomings of SHAM mentioned 
above. We focus on the SFHs of central galaxies, and 
defer a discussion of satellite galaxies to a forthcoming 
paper. 

In general the mass growth/evolution of a central 
galaxy consists of three components: (1) in situ star for- 
mation, (2) accretion of satellite galaxies ('cannibalism'), 
and (3) mass loss due to the evolution of stars. With 
the CMSF model presented in Y12, we can obtain good 
constraints on the growth of the central galaxies, the 
available contribution from satellite galaxies, as well as 
the mass loss due to stellar evolution. With the help of 
the observational constraints on the star formation rates 
(hereafter SFRs) of central galaxies in the local Universe, 
we will be able to obtain the SFHs of central galaxies as 
a function of host halo mass. 

This paper is organized as follows. In Section [2] we 
present the observational data used to constrain the 
SFRs of central galaxies at low redshifts. In Section [3] 
we describe our methodology to constrain the SFRs of 
central galaxies. The results are discussed in Section SI 
where we also present an analytical fitting formula for 
the SFRs of central galaxies as function of halo mass and 
redshift. Model predictions, including the SFHs, stellar 
mass densities, star formation efficiencies, and the frac- 
tion of stars formed in situ, all as functions of redshift, 
halo mass and stellar mass, are presented in Section [SJ 
Finally, in Section [HI we summarize the main findings of 
this paper. 

Throughout this paper, we use the ACDM cosmol- 
ogy whose parameters are consistent with the seventh- 
year data release of the WMAP mission: Om = 0.275, 
Vtjy = 0.725, h = 0.702, and erg = 0.816, where the re- 
duced Hubble constant, h, is defined through the Hub- 
ble constant as = lOO/i km s~i Mpc'^ (WMAP7, 
Komatsu et al. 2011). We use 'In' and 'log' to denote 
the natural and 10-based logarithms, respectively. Un- 
less specified otherwise, throughout this paper we use 
the following units: SFRs are in MQyr~^, SSFRs are in 
yr~^, stellar masses are in /i^^M© and have been derived 
assuming a Kroupa (2001) IMF, and halo masses are in 

2. LOCAL OBSERVATIONAL CONSTRAINTS 

In order to model the SFHs of central galaxies in 
dark matter halos of different masses, we need a reliable 
method to identify central galaxies, as well as accurate 
estimates of their SFRs, preferably across cosmic time. 
At low redshift, such information is available from galaxy 
group catalogs that can be used to represent dark matter 
halos. Here we make use of the SDSS group catalog con- 
structed by Yang et al. (2007; hereafter Y07) using the 
adaptive halo-based group finder developed by Yang et 
al. (2005). The original catalog, based on SDSS Data Re- 
lease 4 (DR4) is updated for DRTQ. The galaxies used for 
the group identifications are selected from the New York 
University Value-Added Galaxy catalog (NYU-VAGC; 
Blanton et al. 2005b), which is based on SDSS DR7 
(Abazajian et al. 2009) but contains a number of im- 
provements in data reduction. From the NYU-VAGC, 

The data is available at |http : //gax . shao . ac . cn/data/Group . html | 





log [M, / h-'Mg] log [M, / h-'M^] 

Fig. 1. — The star formation rates (SFR; right-hand panels) and specific star formation rates (SSFR; left-hand panels) of central galaxies 
as function of their stellar mass A/» (upper panels) and halo mass A/|j (lower panels). The contours in each panel show the number density 
distribution of the galaxies in logarithmic scale, with two neighboring levels differing by a factor two. The shaded vertical histograms in 
each panel are the SFR/SSFR distributions in logarithmic bins of stellar mass or halo mass with widths of ±0.2 dex. The stars in each 
panel indicate the medians of these distributions, while the solid lines in the right-hand panels show our model prediction (see text). 



we select all galaxies in the Main Galaxy Sample with an 
extinction-corrected apparent magnitude brighter than 
r = 17.72, with redshifts in the range 0.01 < z < 0.20, 
and with a redshift completeness Cz > 0.7. The resulting 
galaxy catalog contains a total of 639, 359 galaxies, with 
a sky coverage of 7748 square degrees. A small fraction 
of galaxies have redshifts taken from the Korea Institute 
for Advanced Study (KIAS) Value-Added Galaxy Cata- 
log (VAGC) (e.g. Choi et al. 2010). A total of 36,759 
galaxies, which do not have redshift measurements due to 
fiber collisions, are assigned the redshifts of their near- 
est neighbors. In the present paper, we use the group 
catalog 'modelC, which is constructed on the basis of 
all the galaxies (including those with assigned redshifts) 
and uses model magnitudes to estimate galaxy luminosi- 
ties. In total our catalog contains 472,416 groups, of 
which about 23, 700 have three member galaxies or more. 
Following Y07, we assign a halo mass to each group ac- 



cording to the ranking of its characteristic stellar mass, 
defined as the total stellar mass of all group members 
with ° iMr - 5 log /i < -19.5, where °-'^Mr - 5 log /i is 
the absolute r-magnitude, K- and E-corrected to z = 0.1, 
the typical redshift of galaxies in the SDSS redshift sam- 
ple. The halo mass function adopted in the ranking is 
the model of Tinker et al. (2008), which assumes the 
WMAP7 cosmology and uses A = 200, with A the av- 
erage mass density contrast within the halo (assumed to 
be spherically symmetric). Using the group catalogue, 
we divide galaxies into centrals and satellites; the most 
massive group member is identified as the central galaxy, 
while all other group members are assigned the status of 
satellite galaxy. As already mentioned above, in this pa- 
per we focus on the SFHs of central galaxies only. 

The SFRs and the specific star formation rates, SS- 
FRs (defined as the SFR divided by the stellar mass 
M*), for individual galaxies adopted here are obtained 
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from the data release of Brinchmann et al. (2004 j3, who 
estimated the SFR and stellar mass for each galaxy by 
fitting the observed SDSS spectrum with a spectral syn- 
thesis model. About 94% of the galaxies used in our 
group identifications have estimated SFRs and SSFRs in 
the Brinchmann et al. data release (the vast majority 
of those lacking estimates are fiber collided galaxies with 
assigned redshifts). 

Fig. [T] shows the SSFRs and SFRs of central galaxies 
versus stellar mass M, (upper panels) and halo mass 
(lower panels) . In order to have a volume-limited sample 
of galaxies in stellar mass, we adopt, for given rcdshift z, 
the following stellar mass completeness limit. 



log[A/,ai,„/( /i^^Mq)] = 
4.852-)- 2.2461og£'L(2) 



(1) 



1.123 log(l + z) - 1.186z 



1 - 0.0672 

(see van den Bosch et al. 2008). In addition, we also 
apply the following halo mass completeness limit, 

log[Mft,iim/( h-^Me)] = (z - 0.085)/0.069 + 12 (2) 

(see Yang ct al. 2009b; hereafter Y09b). The contours 
shown in Fig. [1] are the number density distributions 
of central galaxies in 'volume-limited' samples complete 
in both stellar mass and halo mass according to M^,^im 
and Mh,\im given above, i.e. for a given z we only select 
systems with > M^ji^ and > M^^n^^. 

As is evident from the figure, the distributions in both 
SSFR and SFR are bimodal. For given A/, or Mh, 
the central galaxies appear to be separated into two 
distinctive populations, one with high (S)SFRs (here- 
after the "star forming" population) and the other with 
(S)SFRs that are more than 10 times smaller (hereafter 
the "quenched" population). To see this more clearly, we 
show in each panel, using the vertical shaded histograms, 
the distribution of galaxies within given logarithmic stel- 
lar mass (or halo mass) bins with widths of ±0.2dex. 
The star plotted on each of the histograms in each of 
the panels of Fig. [T] indicates the median value of the 
corresponding distribution. 

Let us first focus on the SSFR - stellar mass rela- 
tion shown in the upper left panel of Fig. [T] The 
SSFRs change significantly with stellar mass. Central 
galaxies with stellar masses > 10" /i-^Mq are predom- 
inantly quenched, while those with M* < 10^'^ /i~^Mo 
arc mostly star forming. Note, however, that if satellite 
galaxies were included, a significant fraction of the galax- 
ies at the low mass-end would belong to the quenched 
population (e.g., van den Bosch et al. 2008; Peng et al. . 
2012; Wetzel et al. 2012a, b). Galaxies with intermediate 
stellar masses show strong bi-modal distributions, with 
the quenched population becoming increasingly more im- 
portant as the stellar mass increases. On average the 
SSFRs of central galaxies decrease with increasing stel- 
lar mass. Contrary to the SSFR, however, the total SFR 
distribution depicted in the upper right panel of Fig. [T] 
shows that the SFRs increase roughly linearly with stellar 
mass for the star-forming population, and with a some- 
what slower rate for the quenched population. 

Next, let us look at how SSFR and SFR depend on halo 
mass. As shown in the lower two panels of Fig. [l] galax- 

® see |http : //www ■ mpa-garching . mpg . de/SDSS/DR7/ 1 



ics in halos with masses > 10^^ ° H'^Mq are dominated 
by the quenched population. However, this mode is much 
more 'stretched' in halo mass than in stellar mass. This 
is expected, because the stellar mass of central galaxies 
in massive halos increases only slowly with halo mass, 
M* cx as shown in Yang et al. (2008). 

Finally, we emphasize once more that the above results 
are for central galaxies only. Had we included satellite 
galaxies, the quenched population would have been sig- 
nificantly more prevalent. In what follows these results 
will be used as local observational constraints to model 
the SFHs of central galaxies. 

3. FROM CONDITIONAL STELLAR MASS 
FUNCTION TO STAR FORMATION HISTORIES 

In order to model the SFHs of central galaxies as func- 
tion of halo mass, we first model how the stellar com- 
ponents in dark matter halos evolve with redshift. In 
this section, we first describe our model for the stellar 
mass-to-halo mass relation across cosmic time, based on 
the results published in Y12, followed by a description of 
how such a model can be extended to describe the SFHs 
of central galaxies as functions of halo mass. 

3.1. The conditional stellar mass function of galaxies 
and its evolution 

For a given dark matter halo, the total stellar mass 
it contains can be divided into three components: that 
in the central galaxy, that in the satellite galaxies, and 
that in the form of diffuse halo stars. In Y12 we have de- 
veloped a self-consistent model for the conditional stel- 
lar mass function (CSMF), which describes the stellar 
mass distribution of galaxies in halos of a given mass, 
and its redshift evolution. This model properly takes 
into account that (i) subhalos are accreted at different 
times, and (ii) the properties of satellite galaxies may 
evolve after accretion. Since satellite galaxies were them- 
selves centrals before accretion, the satellite population 
observed today serves as a "fossil record" for the forma- 
tion and evolution of galaxies in dark matter halos over 
the entire cosmic history. Using the observed galaxy stel- 
lar mass functions out to z 4, the conditional stellar 
mass function (hereafter CSMF) at z ^ 0.1, obtained 
from the SDSS galaxy group catalogue, and the two- 
point correlation function (2PCF) of galaxies at z ~ 0.1 
as function of stellar mass, Y12 obtained the relationship 
between galaxy mass and halo mass over the entire cos- 
mic history from z^4toz = 0.1. This relation was then 
used to predict the assembly histories of different stellar 
mass components (centrals, satellites, and halo stars) as 
a function of halo mass. For completeness, we start with 
a brief description of the Y12 results that are relevant 
for the discussion that follows. 

According to Y12, the CSMF of galaxies in halos of a 
given mass can be described by the functional form, 

$(M,|M, z) = $c(Af*|A/, z) + $,(to,|A'/, z) , (3) 

where $c(Af*|A'/, z) and $s("t-*|A/, z) are the contribu- 
tions from the central and satellite galaxies, respectively. 
The CSMF of central galaxies is given by a lognormal 
distribution, 



$c(A/*|A/,z) 



1 



27r(T, 



-exp 



(logAfjAf,,,)^ 



2a2 



(4) 
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where log M,.c is the expectation value of the (10-based) 
logarithm of the stellar mass of the central galaxy and Oc 
is the dispersion (see Y09b). For simplicity, is assumed 
to be independent of halo mass, which has observational 
support (More et al. 2009). Following Y09b, the median 
stellar mass of the central galaxy is assumed to be a 
broken power-law. 



a+/3 



(l + M/./Mi)'^ 



(5) 



so that Af^^c oc M^+'^ (M*,c oc M^) for < Afi 
[Js/'lh ^ Ml). This model contains four free parameters: 
an amplitude Af*_0: a characteristic halo mass. Mi, and 
two power indices , a and /3. All four parameters may 
depend on redshift, as described below. 

The CSMF for satellite galaxies can formally be writ- 
ten as. 



$s(m*|M,„z)^ y dm^y J dMaJ drj 

z 

$e(TO*|ma, Za, z) nsnhirUa, Za\Mh, z) 

PiMa, Za\Mh, z) P{ti) eipt t^i ~ At) , 



(6) 



where 8(a;) is the Heaviside step function and At is 
the time between rcdshifts Za and z. This model 
contains the following ingredients: (i) the accre- 
tion and mass distribution of subhalos, specified by 
nsuhirria, Za\Mh, z) dmadhi(l + Zg) which is defined as 
the number of subhalos in a host halo of mass M^ identi- 
fied at redshift z as a function of their accretion masses, 
ttIq and accretion rcdshifts, Zq (Yang et al. 2011); (ii) 
the growth of the main branch of the host halo, specified 
by P{Ma, Za\Mh, z), which describes the probability that 
the main progenitor of a host halo of mass Mh at redshift 
z has a mass Ma at Za (Zhao et al. 2009); (iii) the or- 
bital distribution of accreted subhalos, specified by P{ri), 
where rj is the orbital circularity (Zenter et al. 2005); (iv) 
the disruption of subhalos during their evolution in the 
host, which is specified by the dynamical friction time 
tdf (Boylan-Kolchin et al. 2008); (v) the relative disrup- 
tion rate of satellite galaxies with respect to subhalos, 
which is characterized by a free parameter pt: pt = oo 
corresponds to no tidal stripping/disruption of satellite 
galaxies, while pt = corresponds to instantaneous dis- 
ruption of all satellites, and finally (vi) the evolution of 
a satellite from its time of accretion until the time cor- 
responding to redshift z, specified by $e('TT-* |ma, Za, z), 
the probability for a subhalo of mass nia and accretion 
redshift Za to host a satellite galaxy of stellar mass m* 
at redshift z. In Y12 it is assumed that 



^cim'Jnia, Za, z) 



1 



27rCT; 



-exp 



(log™*/™ 
2a'2 



(7) 



The dispersion, a'^, is assumed to be the same as (Tc, while 
the median is written as 



(1 - c)m.tM 



cm 



(8) 



where m*_a is the stellar mass of the satellite at the 
accretion time Za, to*.^ is the stellar mass of the cen- 
tral galaxy of a halo of mass nia at time z, and c is 
a free parameter. Thus, if c = then fn'^ = m* a so 



that the stellar mass of a satellite is equal to its orig- 
inal mass at accretion. Physically this corresponds to 
a picture in which a satellite galaxy is instantaneously 
quenched upon accretion. On the other hand, if c = 1 
then m'^ = to*.^ so that the stellar mass of a satellite is 
the same as that of a central galaxy in a halo of mass nia 
at redshift z. This corresponds to the assumption made 
in SHAM; in other words, by setting c = 1 we are ef- 
fectively mimicking SHAM. Note that for c = we have 
that $c(™'* l^a, z) has the same form as the CSMF 
of central galaxies at the accretion redshift Za- 

At any particular redshift, the above model of the 
CSMF is fully described by the following seven free pa- 
rameters: M,^Oj Ml, a, j3, (Tc, Pt and c. In order to 
describe the evolution of the galaxy distribution over cos- 
mic time, Y12 assumed the following redshift dependence 
for the model parameters: 

log[M,,o(z)] =log(M,,o) + iiz 
log[Mi(z)]=log(Mi) +722: 
a{z) = a + 73Z 



log[, 



i[log(/3)+74Z + 752',2] (9) 



(7^(2) = max[0. 173, 0.2z] 
Pt{z)=pt 
c(z) = c . 

The model is thus specified by a total of 11 free param- 
eters, four that describe the CSMF at z = 0, five that 
describe their evolution with redshift, and two (c and pt) 
that describe the evolution of satellite galaxies. 

Once the free parameters are given, the formalism de- 
scribed above allows one to predict the stellar mass func- 
tions as well as correlation functions of galaxies. Thus, 
one can use observational data on the stellar mass func- 
tions (SMFs) and 2-point correlation functions (2PCFs) 
to constrain the model parameters. Y12 used a Markov 
Chain Monte Carlo (hereafter MCMC) method to ex- 
plore the likelihood function in the 11-dimensional pa- 
rameter space, adopting the WMAP7 cosmology. In par- 
ticular, two types of analysis were carried out for com- 
parison. The first uses the SMFs at different rcdshifts 
together with the 2PCFs of galaxies at low-z as con- 
straints; the second uses the SMFs at different rcdshifts 
together with the CSMFs at low z. Since both analyses 
gave very similar results, in this paper we only use the 
results based on the latter analysis. Furthermore, Y12 
used two sets of high redshift stellar mass functions, one 
obtained by Perez-Gonzalez et al. (2008; referred to as 
SMFl), and the other obtained by Drory et al. (2005; 
referred to as SMF2). Since these two data sets reveal 
quite large discrepancies with respect to each other and 
dominate over all the uncertainties in our model ingredi- 
ents, we will present results separately for both of them. 

3.2. The growth of stellar components in dark matter 

halos 

Once the redshift-dependent CSMFs are obtained, one 
can predict the growth of the stellar masses of both cen- 
tral and satellite galaxies along the main branch of their 
dark matter halos. The median stellar mass at redshift 
z of a central galaxy that at redshift zq < z is located in 
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a host halo of mass Mq can be written as 
M,,,(z|Afo,zo) = A4,c(A4,z), 



(10) 



dA^.ioss(z) 
dt 



(13) 



Here, as before, Ma is the median mass at redshift z > zq 
of the main progenitor of a host halo of mass Afg at 
redshift zq. The median totoZ stellar mass of the surviving 
satellite galaxies in the main branch can be obtained by 
integrating the CSMF of satellite galaxies: 

m^,s{z\Mo,zo) ^ y dlog?«* $s(?«*|A/a, z) . (11) 

Thus, once the assembly history of a dark matter halo 
is known, it is straightforward to use the CSMF to ob- 
tain the corresponding assembly histories of the stellar 
components. In addition, one can also estimate the total 
mass of all satellite galaxies that have been accreted into 
the main branch, which includes stellar mass in both the 
currently surviving satellites and those that have been 
cannibalized by the central galaxy or disrupted by the 
tidal field. This is given by 

m^^s,cc{z\Mo,zo) = [ /dlogmH.m* (12) 

$c("T.*|ma, Zo, z) nsy,h{ma,Za\M(), Zq) . 

Note that $0 accounts for evolution in the stellar masses 
of satellite galaxies after accretion due to star formation, 
stellar mass stripping and mass loss due to stellar evolu- 
tion. 

The difference between m,.acc and m,^^ gives the total 
mass of (destroyed) satellites that are either cannibalized 
by the central or disrupted by the tidal field. As shown 
in Y12, the mass in this stellar component in a massive 
cluster can be much larger than that of the central galaxy. 
Hence, a significant fraction of the total stellar mass from 
the disrupted satellite galaxies cannot be associated with 
the central galaxy, but instead has to be in the form of 
diffuse halo stars. Unfortunately, what fraction of stars 
in a halo is associated with such a diffuse halo component 
(also called 'intra-clustcr light' in the case of clusters) is 
still poorly constrained observationally (e.g. Gonzalez et 
al. 2005 Seigar et al. 2007), so that it remains unclear 
what fraction of the 'destroyed' satellite is cannibalized 
by the central vs. added to the stellar halo (e.g. Purcell 
et al. 2007; Kang & van den Bosch 2008; Yang et al. 
2009a). Two extreme assumptions can be made for the 
contribution of the destroyed satellites to the growth of 
the central: (i) minimum (zero) contribution; and (ii) 
maximum contribution, where the contributed mass is 
equal to either the mass growth of the central or the 
total mass of destroyed satellites, whichever is smaller, 
in a given time interval. 

In general, the growth (evolution) of the central galax- 
ies consists of the following three contributions: (i) its in 
situ star formation; (ii) the accretion of stars from satel- 
lite galaxies; and (iii) its passive evolution (mass loss). 
The model described above for the stellar mass assembly 
histories of central galaxies and the possible contribution 
from accreted satellites can therefore be used to infer the 
SFH of central galaxies: 



where 



dA'/,,ioss(z) _ /■*(^UAf,,e(z) d/passivc(t(z) -ii) 



dt 



dti 



At 



dti, 
(14) 

is the stellar mass loss due to the passive evolution 
of stars. Corresponding to the two extreme cases for 
the contribution of the destroyed satellites are two ex- 
treme estimates for the SFH: (i) maximum star forma- 
tion (hereafter 'MAX'), in which the stellar mass growth 
of central galaxies is entirely due to in situ star forma- 
tion, i.e. /acc = 0; (ii) minimum star formation (hereafter 
'MIN'), in which the contribution from accreted satellites 
is maximized (/acc = 1) or SFR{z) = if /acc = 1 leads 
to SFR{z) < 0. 

The stellar mass loss due to passive evolution is ac- 
counted for via the function /passive (Oi which describes 
the mass fraction of stars that at time t after their for- 
mation is still in the form of stars. We obtain /passive (i) 
from the stellar population model of Bruzual & Char- 
lot (2003), kindly provided by Stephane Chariot (private 
communication). Throughout this paper we adopt the 
results for the Kroupa (2001) IMF, which asymptotes to 
/passive ^ 0.57 at latc stages of evolution. For Chabrier 
(2003) and Salpeter (1955) IMFs, these asymptotical 
values are ~ 0.54 and ^ 0.70, respectively. We have 
tested that changing the IMF from Kroupa to Chabrier 
or Salpeter results in changes in the SFRs of < 0.1 dex, 
much smaller than the uncertainties from other sources, 
e.g., the use of SMFl or SMF2. 

4. THE STAR FORMATION HISTORIES OF 
CENTRAL GALAXIES 



SFR(z) . 



dt 



/a 



d[mH.,acc(^) - m^^siz)] 
dt 



Applying the model outlined in Section 13.21 from high 
to low redshifts, we obtain the SFHs of central galaxies 
in halos of different masses. Figs. [2] and [3] show the me- 
dian SFRs and SSFRs as a function of redshift for the 
centrals in halos of different final masses, as indicated in 
each panel. In each panel, the black and blue lines cor- 
respond to the results obtained using SMFl and SMF2, 
respectively. For each set of observational data (SMFl 
or SMF2), predictions based on both the MIN and MAX 
assumptions are shown, connected by the correspond- 
ing shaded As described in Section [321 MIN and 
MAX correspond to the maximum and minimum (zero) 
contribution of satellites to the mass growth of central 
galaxies, respectively, so that they represent the lower 
and upper limits on the SFRs of central galaxies. As one 
can sec, dX z > 2, these upper and lower limits are always 
very similar for halos of all masses. This indicates that 
the growth of stellar mass at these high redshifts is com- 
pletely dominated by in situ star formation (see Y12). 
At lower redshift, the MIN and MAX assumptions re- 
sults in fairly different SFRs, especially for halos with 
Adh > 10^^ /i~^M0. In these halos, mass growth due to 
the accretion of satellite galaxies could be responsible for 
a significant fraction of the final stellar mass of centrals. 
In the absence of accurate estimates for the mass com- 
ponents in halo stars, it will be difficult to discriminate 
between these MIN and MAX models, and therefore to 
tighten the constraints. 
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Fig. 2. — The star formation rates (SFRs) of central galaxies as a function of rcdshift in halos of different present-day masses, as indicated 
in each panel. Black and blue shaded regions reflect the model predictions obtained using SMFl and SMF2, respectively. For each of these 
two cases, the shaded areas mark the SFRs between the MIN and MAX assumptions (see text for details), which become more and more 
similar at higher rcdshift and for lower mass halos. In the panels for halos with Mj, > 10^^ /i^^Mq, local observational constraints from 
SDSS arc shown as the vertical shaded histograms (distributions) and stars (median). The dotted, dashed, and long-dashed lines are the 
MIN, MAX and OBS fits to the SFHs discussed in the text. 
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Fig. 3. — Same as Fig, [2] except that he re we show the specific star formation rates (SSFRs) of central galaxies as function of rcdshift. 
The dashed line is the model fit of Eq. I I15I I. which accurately describes the SSFRs of central galaxies prior to their quenching (see text for 
detailed discussion). 



However, as shown in Section [5J the (S)SFRs of low-z 
central galaxies can be estimated for halos of different 
masses directly from observational data. This provides 
a direct constraint on the SFH at z 0. In Figs. [2] and 
[3] the vertical shaded histograms show the local observa- 
tional constraints for halos with masses > 10^^ /i~^M0. 
Results for halos with lower masses are not available from 
the group catalog (see Y07). In general, the median 
(S)SFRs obtained directly from the SDSS group cata- 
log (indicated by crosses) nicely falls within the shaded 



area, between the upper and lower limits corresponding 
to the MIN and MAX assumptions, indicating that the 
SFRs obtained from our model are consistent with the 
results obtained directly from the SDSS groups. How- 
ever, there are two exceptions: our model under-predicts 
the median (S)SFR of centrals in massive halos with 
Mh > lOi"* " h-^MQ in the case of SMF2, and over- 
predicts the median (S)SFRs of central galaxies in halos 
with Mh = 10^2 /i-iMg. In the case of SMF2, the model 
prediction of the (S)SFR for central galaxies in massive 
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Fig. 4. — Left-hand panel: the peak redshift of the SFR as a function of present-day host halo mass. Data points with error bars are 
the rough estimations of the peak redshifts from the SFHs shown in Fig. [2] with the err orb ars reflecting the redshift interval over which 
the SFIl is within 10% of its peak value. The solid line is a fit to the data points (Eq. |16j V As an illustration, we also show using the 
shaded region where our SFH results are obtained from extrapolations. Right-hand panel: the star formation histories of central galaxies, 
normalized by the stellar masses of the central galaxies at redshift z = 0.1. Over about 4 orders of magnitude in halo mass and stellar 
mass, these normalized SFHs peak at 10~^'^/i^yr~^ (within ~ 0.2dex), as indicated by the horizontal band. 



lialos with Mh > 10^^ ° H^^Mq is well below the median 
(S)SFRs obtained from the group catalog. This happens 
because the use of SMF2 slightly over-predicts the num- 
ber density of very massive galaxies at redshift z ^ 1 
in comparison to that obtained directly from the SDSS 
data. 

For halos with Mh ~ 10^^ /i"^Mq the models based 
on both SMFl and SMF2 predict a median SSFR at 
z = 0.1 of - 10"i"yr-i. In contrast, the SSFRs ob- 
tained directly from the SDSS data reveals two distinc- 
tive populations: a star-forming population, whose SS- 
FRs peak at ~ 10~^°'^yr^^ and a quenched population 
with (SSFR) ~ 10""-^ yr-^ Thus, it appears that our 
model fails to account for the quenched population. We 
suspect that this is caused by two effects. First, due to 
contamination in the group finder, some of the red cen- 
trals in the group catalog arc actually satellite galaxies. 
However, as shown in Y07, this contamination fraction is 
about < 10%, which is not sufficient to explain the rela- 
tively large fraction of quenched centrals in groups with 
Mfi ~ 10^^ /i^^Mq. Second, some of the quenched cen- 
trals may be a population of galaxies that is located in 
the outskirts of a more massive halos which it traversed 
in the past. Star formation in such galaxies is likely to 
have experienced the same quenching as satellites during 
their journeys through the massive halos (Wetzel et al. 
2013). Numerical simulations suggest that such a popu- 
lation is expected, as some of the low-mass halos located 
near massive ones have indeed passed through massive 
halos at least once in the past (e.g. Lin et al. 2003; 
GiU, Kncbc & Gibson 2005; Wang et al. 2009a; Lud- 
low et al. 2009). Observationally, Wang et al. (2009b) 
found that many of the red dwarf galaxies not contained 
within the virial radii of massive halos are located within 
three virial radii from nearby massive halos, consistent 
with the idea that they were quenched when they passed 
through their massive neighbor (see also Geha et al. 
2012). If such a population also exist for galaxies with 



stellar masses ^ 10^" H^'^Mq, it might be possible to 
explain the quenched population of centrals in galaxy 
groups with masses ~ 10^^ h~^MQ. Such a population 
is missed in our model, where all halos that have once 
been accreted into a host halo are considered as subha- 
los hosting satellite galaxies. Clearly, detailed analysis is 
needed in order to quantify the contribution of this kind 
of "passing-through" galaxies to the quenched popula- 
tion (see e.g., Wetzel et al. 2013). 

4.1. Analytical Model 

Let us now move on to a quantification of the predicted 
SFHs. We first focus on the SSFRs of central galaxies. 
As one can see from Fig. [31 at high z the SSFRs roughly 
reveal a power law dependence on redshift. The results 
obtained from SMFl and SMF2 are remarkably similar, 
and are well described by 

log(SSFR/yr^i) = 2.51og(l-hz) (15) 
-0.121og(M,,//i-iMo) - 12.0, 

which is indicated as the dashed line in each of the panels 
in the figure. In smah halos with Mh < 10^^ ° /i^^Mq, 
Eq. ([151) liolds for the entire SFH. In halos with Mh > 
2^Qi2.5 /j-iMq, however, Eq (UHl) holds only for z > Zpk, 
where Zpk is defined as the redshift at which the SFR 
peaks (see below). Eq. (|15p indicates that the SSFR of 
central galaxies depends only weakly on galaxy mass. 
For example, central galaxies in halos with Alh = 10^^ 
and 10^2 h-'^MQ have similar SSFRs, even though their 
stellar masses differ by more than an order of magnitude 
(cf. Fig. 13 in Y12). The dependence of SSFR on red- 
shift is much stronger, increasing by a factor of about 6 
from z = 1 to z = 3. Interestingly, this scaling of the 
SSFR with halo mass and redshift is almost identical to 
that of the specific mass accretion rate of dark matter 
halos, which scales as Mh/Mh oc M^-^^ (1 + z)^-^^ (Dekel 
et al. 2009; see also McBride et al. 2009; Fakhouri et 
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al. 2010 for similar results obtained from N-body sim- 
ulations). This suggests that the SFR of star-forming 
(i.e., non-quenched) central galaxies is regulated by the 
rate at which their host halos accrete mass, in excellent 
agreement with a number of recent studies (e.g.. Button, 
van den Bosch & Dekcl 2010; Bouchc ct al. 2010; Dave, 
Finlator & Oppenheimer 2012). 

The SFRs in Fig. [2] show that Zpk increases with in- 
creasing stellar mass (or halo mass), indicating that more 
massive centrals are quenched earlier (a manifestation of 
what is often called 'downsizing'). We have estimated 
Zp\^ as a function of halo mass, Mh, using the SFHs ob- 
tained from both SMFl and SMF2. The open squares in 
the left-hand panel of Fig. |4] show the average between 
SMFl and SMF2, while the errorbars arc obtained from 
the points where the SFH is about 10% of the peak value. 
Using the following functional form to fit the data, 

Zpk = max[a(log Mh - 5), 0] , (16) 

we obtain a = 0.568 and b — 10.10. This fitting function 
is shown in the left-hand panel of Fig. U] as the solid line. 

In addition to Zpk increasing with halo mass, the actual 
SFR at the peak redshift also increases with halo mass. 
As is evident from Fig. [51 this halo mass dependence be- 
comes weaker in more massive halos. This kind of halo 
mass dependence is also seen for the luminosities and 
stellar masses of central galaxies at low redshift (see, e.g., 
Y09b). Interestingly, if we normalize the SFRs of the cen- 
tral galaxies by their stellar masses at 2 = 0.1, the peak 
amplitude becomes virtually independent of halo mass 
and stellar mass! This is demonstrated in the right-hand 
panel of Fig. HI which plots log[S'Fi?/M*.o.i] as func- 
tion of redshift for five different values of the present-day 
halo mass, as indicated. The peak amplitude is about 
log[SFR/M*_o.i X /i~^yr] = —9.3, which is indicated by 
the horizontal line. The hatched, horizontal band cov- 
ers ±0.2 dex around this line and roughly represents the 
uncertainty in the estimate of the SFR and stellar mass. 
Thus, central galaxies with stellar mass M^.o.i at z = 0.1, 
have a maximum SFR 

SFRpk _ M*,o.i , s 

[Mgyr-i] " 10«-3/i-2]vig ^ ) 

at redshift Zpk given by Eq. p6)) . These equations hold 
for central galaxies covering over 4 orders of magnitude in 
both halo mass and stellar mass. Note that our modeling 
here is based mainly on results obtained from SMFl. The 
results obtained from SMF2 are qualitatively similar, and 
shown in the Appendix for comparison. 

Motivated by the general appearance of the SFHs, we 
use the following functional form to model the SFHs of 
central galaxies as function of halo mass and redshift: 

SFR(M., .) = SFRpkXexp {- "''^^'l:^^^^^ '^'^^ } , 

(18) 

where ^(zpk) describes the decay of the SFR with respect 
to the peak. By trial-and-error, we find that a simple 
power- law scaling with (1-1- Zpk) can adequately describe 
the SFHs. For z > Zpk we obtain 

a(zpk) = 0.0576(1 + 2pk)°-™' , (19) 
where the parameters are obtained using least squares 



fitting to the SFHs at high rcdshifts. The differences 
between the SMFl and SMF2 model predictions are used 
as weights in the fitting. For z < Zpk we obtain different 
power-law scalings: 

( 0.0762(1 + Zpk)°-^^3 (OBS) 
a(zpk) = < 0.0706(1 + Zpk)°-^*° (MAX) . (20) 
[ 0.317(1 + Zpk)-2-i° (MIN) 

Here the case marked 'OBS' assumes that the SFHs at 
low-z are constrained by the SFR measurements from 
the SDSS groups using least squares fitting to the me- 
dian SFRs of central galaxies as a function of stellar mass 
(asterisks in upper right-hand panel of Fig. [1]). The case 
marked 'MAX' assumes that all stars in central galax- 
ies are formed in situ, in which case the stellar mass of 
central galaxies as function of redshift, zo, follows from 
integrating the SFH and taking into account the effect 
due to passive evolution, 

M^zo) = f " SFRMAx(i) /passivo(i.o -t)dt. (21) 
Jo 

Here tz^ is the age of the universe at redshift zq and 
1 — /passive (Af) is the mass fraction of stars formed that a 
time At later has been returned to the IGM due to stel- 
lar (passive) evolution. By definition, M*(zo = 0.1) = 
-^^*,o.ij which is used to obtain the two parameters in 
Eq.' ipU)) for the 'MAX' case. For the 'MIN' case, the two 
parameters in Eq. (|20p arc obtained from the correspond- 
ing low-redshift SFHs given by SMFl. 

The resulting SFH model predictions are plotted as 
the long-dashcd (OBS), dashed (MAX) and dotted lines 
(MIN) in Fig. O Thus, Eq. ^ describes the median 
SFH for central galaxies with a given stellar mass at red- 
shift z = 0.1. For halos of a given mass Mh, one can first 
obtain the corresponding M*_o.i using Eqs. ^ and 
and then obtain the SFH from Eq. ([T5)l . The uncertain- 
ties in the SFH model may be gauged from the differences 
between the model predictions for SMFl and SMF2, and 
between the results for the two extreme assumptions, 
MIN and MAX (see Fig. [2]). As pointed out in Y12, 
at the present the uncertainties in the CSMF modeling 
are dominated by systematic errors between SMFl and 
SMF2 (see Fig. 14 in Y12), rather than by the statisti- 
cal errors in the data. Note that our model is obtained 
from observational measurements of stellar mass func- 
tions that are limited in both redshift (z < 4 for SMFl 
and z < 5 for SMF2) and stellar mass (see Fig.3 in Y12). 
Results beyond these redshift and/or stellar mass ranges 
are in general obtained by extrapolation, and therefore 
less reliable. As an illustration, we show, in the left panel 
of Fig. [4] using the shaded region, where our SFH results 
are obtained from extrapolations. 

Finally, we emphasize that our model for the SFHs 
is based on the assembly histories of central galaxies in 
halos of a given final mass, and it is expected to work well 
only for centrals with mass assembly histories close to 
the median. Deviations from the median are expected to 
result in deviations of the inferred SFHs from the model 
prediction and may in fact be the main source of the 
scatter seen in Fig. [T] In a forthcoming paper, we will 
come back to the modeling of this scatter together with 
the SFHs of satellite galaxies. 
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Fig. 5. — The SFHs of central galaxies, similar to Fig.[2]but here our model predictions are compared with the results obtained in a few 
recent papers. The lower and upper boun dar ies of the shade areas are our OBS model pr edic tions of the SFHs for central galaxies based 
on SMFl with the SFR described by Eq. ^8^, and SMF2 with the SFR described by Eq. [KM . respectively. The dashed, dotted and long 
dashed lines in each panel are model predictions based on Leitner & Kravtsov (2011), Leitner (2012) and Moster et al. (2013), respectively. 
See text for a detailed discussion. 



4.2. How to Use the Model 

As a summary, the following is a brief description of 
the procedure for obtaining the median SFH for a central 
galaxy in a halo of a given mass: 

1. Start from a halo with mass Mh at low redshift 
(e.g. z - 0.1). Use Eqs. (17) and (40) in Yang 
et al. (2012), which are the relations adopted in 
this paper, or use Eq. (20) of Yang et al. (2009), 
which was obtained directly from the SDSS group 
catalog, or use any other stellar mass - halo mass 
relation for central galaxies, to obtain the stellar 
mass, A/*^o.ii of the central galaxy at z ~ 0.1; 



2. Use Eq. (|T6|) to obtain the redshift, Zpk, at which 
the SFR peaks; 

3. Use Eq. (Ull) to get the peak value of the SFR, 
SFRpk; 

4. Finally, the median SFH of the central galaxy is de- 
scribed by Eq. (fTS]) with (T(zpk) given by Eqs. (fT9|) 
and dini) for the case labelled 'OBS'. 



4.3. Comparison with some recent results 

As pointed out in the introduction, there are a number 
of recent studies that, similar to this paper, modelled the 
SFHs of galaxies as function of halo mass. All these stud- 
ies make use of SHAM (or a modified version thereof), 
together with halo merger trees generated from N-body 
simulations, to make predictions for the mass growth and 
SFH of galaxies (e.g., Moster et al. 2013; Behroozi et al. 
2012, 2013; Wang et al. 2012). In ah these analyses, 
assumptions have to be made regarding the evolution of 
satellite galaxies, in particular regarding the contribu- 
tion of satellite accretion to the mass assembly of central 



galaxies. Although the results from all these analyses 
are qualitatively similar to each other and to our results, 
there are substantial, quantitative differences. In this 
subsection, we make a quantitative comparison between 
our results and those obtained by Moster et al. (2013). 

Fig. [5] shows our model predictions of the SFHs for 
central galaxies. The lower and upper boundaries of the 
shaded areas are our model predictions based on SMFl 
(the same as the long dashed lines in Fig. [2]) and SMF2 
(the same as the long dashed lines in Fig. [T5|) . respec- 
tively. The long-dashed line shown in each panel is the re- 
sult obtained with the fitting formula provided by Moster 
et al. (2013). The Moster et al. results are obtained by 
assuming a WMAP7 cosmology and a Chabrier (2003) 
IMF; these arc sufficiently similar to the cosmology and 
IMF adopted here, that a direct comparison is justified. 
Within the uncertainty (~ 0.2 dex in the SFH), the peaks 
in the SFHs of Moster et al. (2013) are roughly consis- 
tent with our model predictions. However, at high-z their 
SFHs are higher than ours while at low-z and for mas- 
sive halos, the Moster et al. SFHs are lower than our 
model predictions. These differences, especially those at 
low-z, most likely comes from the different treatment of 
satellite galaxies. Moster et al. (2013) assume that 20% 
of the stellar mass of accreted satellite galaxies 'escapes' 
from the central galaxy, ending up as diffuse halo stars. 
If we add such a prescription to our model, we predict a 
\ow-z behavior that is very similar to that of Moster et 
al., although we find that a 'satellite disruption fraction' 
of ^ 45%, rather than 20%, yields results that are in 
better agreement with observational constraints on the 
SFHs. 

It is also interesting to compare our results to those 
of Leitner & Kravtsov (2011) and Leitner (2012), who 
modeled the stellar mass growth in star-forming galaxies 
using the observed SFR-stellar mass relations at differ- 
ent redshifts. The short-dashed and dotted lines in Fig. [5] 
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Fig. 6. — Maps of our model predictions for the star formation rates (SFR; upper panels) and specific star formation rates (SSFR; lower 
panels) of central galaxies as functions of halo mass (left-hand panels) or stellar mass (right-hand panels) and redshift. The vertical bars 
on the right-hand side indicate the corresponding color-coding, where the SFR and SSFR are expressed in units of MQyr"'^ and yr~^, 
respectively. The solid curves in each of the panels show the median growth of halo mass (left-hand panels) and stellar mass (right-hand 
panels) along the main branch of halos with different present-day masses. 



show the model predictions of Leitner & Kravtsov (2011) 
and Leitner (2012), respectively. Although their models 
predict SFH peaks that are very similar to our model 
predictions, their models predict SFRs at low (high) red- 
shifts that are significantly higher (lower). The differ- 
ence at low-z is expected, because their models only de- 
scribe star-forming galaxies, while our models are for the 
general population. At high-z, the results of Leitner & 
Kravtsov (2011) and Leitner (2012) arc not obtained di- 
rectly from observational data but rather from extrapo- 
lations of the low-z data to higher redshifts. We there- 
fore suspect that the difi'erence between their results and 
ours reflects an amplification of the uncertainties in their 
SFR-stellar mass relations at lower z. 

5. MODEL PREDICTIONS 

In this section we use our model for the SFHs of cen- 
tral galaxies as function of halo mass to make a num- 
ber of predictions that provide valuable insight into 



how star formation proceeds in the Universe. Note 
that throughout this section, unless specified otherwise, 
our model predictions are made from observational con- 
strained SFHs obtained from SMFl (as SMFl shows bet- 
ter agreement with the SFR measurements at z = 0.1 
than SMF2). 

5.1. The star formation rate as a function of halo 
mass and redshift 

Let us first look at the SFR maps in the Mh - z and 
- z spaces, shown in the upper row of panels of Fig.[6l For 
reference, the solid lines show the growth of the host halo 
mass (left-hand panel) and the growth of stellar mass 
of the central galaxies (right-hand panel). These rela- 
tions are obtained from the halo mass accretion model 
of Zhao ct al. (2009) and from the median stellar mass - 
halo mass relation for central galaxies obtained by Y12 
(Eq. [To]), respectively. Results are only shown for sys- 
tems that have Mh < 10^^ h~^MQ, since more massive 
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Fig. 7. — Model predictions for the star formation rates of central galaxies. The upper panels show the SFRs of central galaxies as 
function of rcdshift for different halo masses (upper left-hand panel) and different stellar masses (upper right-hand panel), as labelled. The 
lower panels show the SFRs of central galaxies as functions of halo mass (left-hand panel) and stellar mass (right-hand panel) at different 
redshifts, as labelled. All these curves are either vertical or horizontal cuts through the SFR maps shown in the upper panels of Fig. [6] 
In each of the four panels, the vertical cut-offs in the curves correspond to median main branch mass of a halo with present day mass 



halos arc extremely rare. As one can see, the SFRs peak 
in massive galaxies with ^ 10^^ ° H^'^Mq in halos 
with Mh 10"-^ /i-^Mq at rcdshift z - 2.5, and the 
peak star formation rate is about lOOMQyr"^. At low 
rcdshift (z < 0.5), SFRs are highest in galaxies with 
M* ~ 10i° ° /i-^Mq in halos with Mh ~ lO^^ h'^MQ. 
Also there is a lower halo mass limit at < 10^°'^ h~^MQ, 
below which the star formation is low over the entire red- 
shift range. Note that these results are in good qualita- 
tive agreement with the results of Behroozi et al. (2012). 

In addition to the SFR, in the lower panels of Fig. IH] 
we show the corresponding SSFR maps. In general the 
SSFR maps show a smooth gradient in rcdshift, being 
higher at higher rcdshift. At high z, the SSFR depends 
only weakly on halo mass and galaxy mass. Only after 
the peak rcdshift of the SFR (see Eq. [H]), does the 
SSFR decrease more rapidly in more massive galaxies 



hosted by more massive halos. 

To better illustrate the behavior of the SFRs, it is use- 
ful to show some horizontal and vertical cuts of the SFR 
maps. The upper panels of Fig. [7] show the SFR as a 
function of redshift for central galaxies in different mass 
halos (upper-left panel) and for central galaxies of differ- 
ent stellar masses (upper-right panel). In the lower pan- 
els, we plot the SFR as a function of halo mass (lower- left 
panel) and stellar mass (lower-right panel) at different 
redshifts, as indicated. The sharp cut-offs in each panel 
correspond to the median main branch mass of a massive 
halo with Mh = 10^^ H^'^Mq at z = 0, beyond which the 
abundance of systems is negligibly small. 

A number of interesting characteristics are evident 
from these plots. First, for central galaxies in massive 
halos with Mh > 10^^ h~^MQ, the SFR increases mono- 
tonically with redshift (see the upper-left panel of Fig. [7]). 
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Fig. 8. — Maps of our model predictions for the star formation rate density (SFRD) of central galaxies. The color coding used is indicated 
in the vertical bar on the right-hand site. These SFRDs are normalized to be the mass density of stars formed per unit log Mh per unit 
log(l + z) (loft-hand panel) or per unit log M» per unit log(l -|- z) (right-hand panel), and are expressed in units of MQyr~^[ h~^yipc]~^ . 
The contours in the two panels enclose a region that contributes 50% of the total stars formed in the universe. As in Fig. [61 the solid curves 
show the median growth of halo mass (left-hand panel) and stellar mass (right-hand panel) along the main branch of halos with different 
present-day masses. 



In less massive halos, however, the SFR-z relation is 
not monotonia. Rather it reveals two maxima, one at 
high redshift (z > 5), and one at intermediate redshift 
z ~ 1. Contrary to the SFRs in more massive halos, the 
SFR-z relation has a local minimum at a redshift that 
increases with decreasing halo mass. However, we cau- 
sion that such two maxima features in small halos are 
obtained from the extrapolations of the current observa- 
tional stellar mass and redshift limits (see left panel of 
Fig. 2]). Second, for central galaxies of given stellar mass 
the SFRs increase monotonically with increasing redshift 
over the entire redshift ranges probed, with a weak indi- 
cation that the rate of increase is larger for more mas- 
sive centrals (see upper- right panel of Fig. [7]) . Third, at 
z > 2 the SFR is higher for more massive halos (see the 
lower-left panel of Fig. [T]). At z < 2, though, the SFRs 
in massive halos are strongly suppressed relative to those 
at higher z. The mass scale at which this suppression be- 
comes apparent shifts to lower halo masses with decreas- 
ing redshift; for smaU halos with < 10^^-^ /i~^Mq, 
the SFR is almost independent of redshift. Finally, as 
is apparent from the lower right-hand panel of Fig. [71 
the SFRs for galaxies of a given stellar mass generally 
increase with redshift, and, for a given redshift, the SFR 
roughly shows a power-law dependence on stellar mass, 
especially at high redshift. These features are quite ro- 
bust for the observationally constrained SFHs, regardless 
whether they are constrained from SMFl or SMF2 (see 
Appendix A). 



5.2. Star formation rate density 

Using the 'SFR maps' shown in the upper panels of 
Fig. [51 we can obtain similar maps but for the star for- 
mation rate density (SFRD). These have the advantage 
that they highlight when and where the majority of stars 
in the Universe formed. The SFRD is related to the SFRs 



according to 

SFRD(A4, z) = InlO x n{Mh, z) M,, (22) 

><^^^(^^-^)dWTiy 

where n(Mfi,z) = dn(z)/dMfi is the comoving number 
density of halos at z with masses in the range [Mh, Mh + 
dMh]. Thus defined, the SFRD is the mass density of 
stars formed per logM/i per log(l + z), and the result is 
shown in the left-hand panel of Fig. [51 The SFRD peaks 
at redshift z - 1 in halos of Mh ^ lO^^ h-'^Mq. The 
solid, green contour shown on top of the SFRD map en- 
closes the region that contributes 50% to the total SFH 
in the Universe, and shows that half of the total stel- 
lar mass is formed in halos with masses in the range 
10" ! H-^Mq < Mh < 1012-3 /j-IMq and in the red- 
shift range 0.4 < z < 1.9. Overall, halos with masses 
Mh = I0i°-5~i3-^ h~^MQ contribute the vast majority of 
all star formation in the universe. At z < 1.5, the star 
formation density becomes progressively more dominated 
by low-mass halos. 

A SFRD map can also be constructed in the stellar 
mass - redshift plane, 

SFRD(Af„ z) = SFRD(Af;„ z)^^^ , (23) 

d log 

and the results are shown in the right-hand panel of 
Fig. [H Here we see that the SFRD peaks at Af* ~ 
10 h~^MQ, and that half of the stars in the Universe 
is formed in galaxies within a narrow (~ 1.3dex) stellar 
mass range: lO''-^ H-^Mq < AT* < lO^^-^ /i-^Mq. 

5.3. Star formation efficiency and gas mass depletion 

rate 

Next we look at the gas mass depletion rates (GMDR) 
of central galaxies in halos of different masses. Since our 
current modeling does not include any gas components. 
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Fig. 9. — Model predictions for the gas mass depletion rates (GMDR; upper panels) and star formation efficiencies (SFE; lower panels) 
of central galaxies as functions of halo mass (left-hand panels) or stellar mass (right-hand panels) and redshift. The color codings used are 
indicated in the vertical bars on the right-hand site. Whereas the SFEs are dimensionless, the GMDRs arc in units of yr~^. As in Fig. [SI 
the solid, black curves show the median growth of halo mass (left-hand panel) and stellar mass (right-hand panel) along the main branch 
of halos with different present-day masses, while the solid, green lines mark the maximum GMDRs/SFEs as function of redshift. 



we define the GMDR to be the SFR normahzed by the 
total baryonic mass, defined as tlie lialo mass times the 
universal baryon fraction fb'. 

GMDR(A4,.) = 55»i^, (24, 

Jb Mh 

where we adopt fb = ilh/^m — 0.167 (e.g., Komatsu et 
al. 2011). Thus defined, the GMDR is the reciprocal of 
the time it takes a galaxy to consume all the gas associ- 
ated with its halo at the current SFR. The top left and 
right panels of Fig. [9] show the GMDR maps in the halo 
mass vs. redshift and stellar mass vs. redshift spaces, 
respectively. The green solid line shown in each panel 
marks the maximum GMDRs as function of redshift. 
Three trends are worth noting. First, the GMDRs in low 
mass halos {Mh < 10^^ h~^MQ) are always strongly sup- 
pressed with respect to their more massive counterparts. 
Second, at z > 3.5, the most massive halos have the high- 



est GMDRs, which are roughly constant at ~ 10^^'^ yr^^ 
over a large range in redshift. These high GMDRs result 
in rapid growth of the central galaxies, as is evident from 
the black solid lines shown in the right-hand panel. Fi- 
nally, below z ~ 3.5 'downsizing' kicks in, in that the 
peak in the GMDR shifts to lower halo masses with de- 
creasing redshift. There appears to be some 'quenching' 
mechanism (or at least some mechanism that manages to 
strongly suppress the GMDR), which operates in halos 
whose mass shifts down with decreasing redshift; cen- 
trals that end up in more massive halos at z = quench 
at a higher redshift, when their halo mass is more mas- 
sive. Whereas the present-day quenching mass is close 
to 10^^ /i~^Mq, centrals that have ended up in the most 
massive halos quenched around z ^ 3.5 when their main 
progenitor halo had a mass ^ 10^^ h~^MQ. 

An insightful way of expressing the SFHs of central 
galaxies is to compare them with the mass accretion his- 
tories of their halos. To this end, wc follow Bchroozi et 
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Fig. 10. — Model predictions for the star formation efficiencies of 
central galaxies as function of the main branch mass, Ma, of their 
host halo. Results are shown for host halos with present-day masses 
log(M^/h-iM0) = 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, 14.0, 14.5, 15.0. 
The present-day mass for each curve is evident from its end-point 
at the high-mass end. 

al. (2013) and define the star formation efficiency (SFE) 
of central galaxies in the main branch progenitors as 



SFE(M;„z) 



SFR(A/fe,z) 

d(/bMa)/di 



(25) 



Thus defined the SFE describes the SFR in the main 
branch progenitor at redshift z (where the halo mass is 
Mg) in units of the baryonic mass accretion rate at the 
same redshiffl The lower panels of Fig. [9] show maps 
of the SFE in the halo mass vs. redshift (left-panel) and 
stellar mass vs. redshift (right-panel) spaces. The SFE 
peaks in halos with mass ^ lO^"'^ "' /i~^M0 at redshift ^ 
0.5, with central galaxy mass ~ lO^'^/i~^M0. Compared 
with the GMDR maps shown in the upper panels, the 
peak in the SFE shifts to lower halo/stellar mass and 
lower redshift. This is simply a reflection of the fact that 
low mass halos have lower specific mass accretion rates 
than their more massive counterparts (i.e., more massive 
halos assemble later). 

These results are in qualitative agreement with those 
obtained by Behroozi et al. (2013), but quantita- 
tively there are differences. In particular, compared to 
Behroozi et al. (2013), our model seems to predict signif- 
icantly stronger redshift dependence. For example, ac- 
cording to our model the halo mass at which the SFE 
is highest increases from ^ 10^^'^ /i^^M© at z = 
to Af/j ^ 10^^'^ ft.~^MQ at z = 3 (see solid, green line in 
lower left-hand panel of Fig. |9l). In contrast, the corre- 
sponding SFE peak halo masses obtained by Behroozi 
et al. (2013) are 10"-^ /i-^Mq and lOi^-^J at z = 
and z = 3, respectively. Interestingly, we find that the 
stellar mass of the central galaxies for which the SFE 
(or GMDR) is maximal is virtually independent of red- 
shift at M*,c ^ lO^ "^ /i-^Mq, but only for z < 2.5. Our 
model suggests that for z > 2.5, this characteristic stel- 

® Note that we are unable to distinguish whether star formation 
at a given redshift consumes previously or newly accreted gas. 



lar mass rapidly increases to effectively become equal to 
that of the most massive galaxies present at those red- 
shifts. Such a feature is reconfirmed very recently using 
a different approach by Lu et al. (2013 in preparation). 

To highlight some of the features of the SFEs, 
Fig. [ini shows cuts of the SFE maps along some of 
the main branch histories (solid lines in the left pan- 
els of Fig. |9l). Results are shown for central galaxies 
with present-day halo masses of log{Mh/ h~^MQ) = 
11.0,11.5,12.0,12.5,13.0,13.5,14.0,14.5,15.0. As the 
main progenitor mass increases with time, the SFE ini- 
tially increases, reaches a maximum, after which it de- 
creases rapidly to a quenched state. Interestingly, the 
initial increase of the SFE with increasing Ma is much 
steeper for central galaxies that end up in less mas- 
sive halos. A rapid increase of SFE over the range 
10^° - 10" /i"^Mq in Ma, which is predicted by our 
model for centrals that end-up in present-day halos with 
Mil < 10^^'^ h~^MQ, is consistent with the presence of 
a 'halo mass floor' Mmin ^ 10^^ h~^MQ, below which 
star formation is strongly suppressed, as suggested by 
Bouche et al (2010). However, our model predicts that 
central galaxies that end up in more massive halos have 
fairly high SFEs when their main progenitor mass is 
~ 10^° Ii'^Mq. This suggests that this halo mass floor 
must have been substantially lower (or absent) at higher 
redshifts (z > 5). 

Another interesting feature of the SFE(A/a) curves in 
Fig. [To] is the evolution in the 'quenching mass', which 
we define as the main-progenitor mass at which the SFE 
is maximal. Present-day halos with Mh > 10^'^ h~^MQ 
all seem to quench when Ma ^ 10^^'^ h~^MQ, at which 
point their SFE is ~ 0.1 (i.e., the star formation rate is 
ten percent of the baryon accretion rate) . For present- 
day halos with Mh ^ lO^'^ Ii'^Mq our model predicts a 
'downsizing' behavior, in that the quenching mass shifts 
to lower Ma for present-day halos that are less massive. 
In addition, the peak value of the SFE increases, coming 
close to unity for Mh ~ 10"-^ /i-^Mq. 

An important, outstanding question in galaxy for- 
mation is what physical process is responsible for the 
quenching of central galaxies, which seems to happen 
whenever the halo mass is of order 10^^ /i~^Mq. Inter- 
estingly, this mass scale is very similar to the one that 
separates cold-mode accretion and hot-mode accretion 
(Birnboim & Dekel 2003; Keres et al. 2005), suggesting 
that the quenching of star formation in central galax- 
ies may be related to the ability of a dark matter halo 
to form a hot gaseous halo. This is indeed what seems 
to be needed to explain the observed bimodality in the 
distributions of galaxy colors and specific star forma- 
tion rates (e.g., Cattaneo et al. 2006; Birnboim et al. 
2007). In addition to cold mode/hot mode accretion, 
other mechanisms have also been invoked to explain the 
quenching of star formation in massive central galaxies, 
ranging from AGN feedback (e.g.. Tabor & Binney 1993; 
Ciotti & Ostriker 1997; Croton et al. 2006; Bower et 
al. 2006; Hopkins et al. 2006) and gravitational heat- 
ing (e.g., Fabian 2003; Khochfar & Ostriker 2008; Dekel 
& Birnboim 2008; Birnboim & Dekel 2011) to thermal 
conduction (e.g., Kim & Narayan 2003). Although it re- 
mains unclear which of these processes dominates, and 
how exactly they operate, it is clear that any successful 



16 



Yang ct al. 



model has to be able to explain why star formation in 
galaxies is quenched once their halo masses reach a char- 
acteristic mass, which 'downsizes' from ~ 10^^'^ h~^MQ 
at high rcdshifts {z > 3.5) to ~ lO^"'^'^ at the present-day. 

5.4. Comparing SFE with central-to-host halo mass 

ratios 

Since the SFE indicates the fraction of newly accreted 
baryonic matter that is converted into stars, the quantity 
SFEx/f, can be regarded as the instantaneous mass ratio 
between the newly formed stars and newly accreted dark 
matter. It will be interesting to compare this mass ratio 
to that between the stellar mass of a central galaxy and 
the mass of its dark matter host halo. The latter has been 
extensively studied in recent years, and is an integration 
of SFEx/b over cosmic time plus a contribution due to 
the accretion of satellite galaxies (important in massive 
halos at low redshift) and the impact of passive evolution 
(which removes ~ 40 percent of the stars formed at early 
times from the mass budget of present day stars). 

We first show, in Fig. [TTl the c/Mh ratios obtained 
in Y12 for SMFl (red shaded curves) and SMF2 (black 
shaded curves). The shaded areas refiect the 68% con- 
fidence levels resulting from the statistical errors in the 
resulting M^,^c/M^h ratios. The differences between the 
SMFl and SMF2 curves can be regarded as roughly re- 
fiecting the systematical errors in the current data. Com- 
pared to the "no-evolution" model, which simply is the 
M^,^c/Mh a.t z = 0.1 (cyan, dashed curve in each panel), 
the Mtf^c/Mh ratio docs evolve significantly, especially 
beyond redshift z ~ 1. On average, the peak Mt,^c/Mh 
ratio at redshift z > 2.0 is ^ 0.5dex below the one at 
redshift z = 0.1. 

The solid line in each panel of Fig. [11] indicates our 
model prediction for the instantaneous ratio, SFE x /t,, at 
the corresponding redshifts, reduced by 40% in order to 
(roughly) correct for passive evolution (see discussion in 
tj3.2[) . Again, for comparison, we also show in each panel 
the "no-evolution" model predictions (long-dashed green 
curve), which simply is the SFE x ratio at z = 0.1. 
Note how the peak of the SFE x vs. Mh curve first 
increases by about -|-0.4dex from z = 0.1 to z ~ 0.6, 
after which it gradually decreases to about — 0.5dex at 
z > 3.0. 

Comparing the 'integrated' mass ratios, Ms,_c/^ih to 
the 'instantaneous' ratios, SFE x /;,, in different redshift 
bins, one notices that the two ratios are in good agree- 
ment with each other at high redshifts {z > 3). This 
is as expected, since (i) the SFE does not show strong 
time evolution at high redshift, and (ii) the contribution 
from satellite galaxies is negligible. Moving to lower red- 
shifts, wc see an ever increasing 'lag' between the 'inte- 
grated' and 'instantaneous' mass ratios at the high mass 
end. This is a manifestation of the quenching of central 
galaxies in massive haloes; their instantaneous SFRs are 
a poor indicator of their integrated (past) SFHs. In ad- 
dition, massive centrals may have accreted a significant 
fraction of their stellar mass in the form of satellite galax- 
ies (see §5.6l below). which is also not reflected in their in- 
stantaneous SFRs. Finally, at all 0.2 ^ z < 2.5, the peak 
in the 'instantaneous' mass ratio is higher than that in 
the 'integrated' mass ratio. This is a manifestation of the 
fact that the instantaneous SFE has a global peak at red- 
shift z ~ 0.5. Hence, at redshifts above and in a certain 



range blow this peak redshift, the instantaneous mass ra- 
tio should be larger than its time-integrated equivalent. 

5.5. The cosmic star formation densities 

Once we know how the SFRs of central galaxies depend 
on Mh and z, we can combine the dependence with the 
halo mass function to predict the star formation history 
of the universe, as described by the cosmic star formation 
density (SFD) deflned as 

/>oo 

SFD(z)=/ SFR(A/,„ z) n(Af,„ z)dA/,, . (26) 
Jo 

Note that our model prediction only accounts for the 
contribution due to central galaxies, whereas the data on 
the cosmic star formation history includes contributions 
from both centrals and satellites. However, since satellite 
galaxies contribute less than ^ 40% of the total galaxy 
population (e.g., Mandclbaum ct al. 2006; Tinker ct al. 
2007; van den Bosch ct al. 2007; Yang et al. 2008; 
Cacciato et al. 2013) and have significantly lower SFRs 
(e.g., Balogh et al. 2000; Weinmann et al. 2006, 2009; 
Wetzel, Tinker & Conroy 2012a, b), their contribution to 
the SFD is sufficiently small that a comparison of our 
model prediction with data is still meaningful. 

In the left panel of Fig. [121 we compare our model 
prediction of the cosmic SFD to a compilation of data 
taken from Hopkins & Beacom (2006). Hopkins & Bea- 
com (2006) converted all the SFRs to a Salpcter (1955) 
IMF, which is different from the Kroupa (2001) IMF 
used in this work. As suggested in Perez-Gonzalez ct al. 
(2008), the stellar masses based on the Salpetcr (1955) 
IMF are systematically larger by a factor of ^ 1.7 than 
those based on the Kroupa IMF. Here we have applied 
such a conversion in our comparison by multiplying our 
model prediction by a factor 1.7. The solid and long 
dashed lines show our model predictions based on SMFl 
and SMF2, respectively, with the latter obtained from 
the fitting formula provided in Appendix A. For com- 
parison we also show, using a shaded band, the range of 
the SFDs obtained from the three models for the SFR, 
'OBS', 'MAX' and 'MIN', and using SMFl and SMF2. 
This band therefore roughly captures the uncertainties 
in our model. 

As one can see, at low redshift (z < 2.0), our model 
prediction is in good agreement with the data. At 
z > 2.0, however, our model under-predicts the cosmic 
SFD compared to the data. Several sources might con- 
tribute to this discrepancy. First, our model prediction of 
the SFD is based on the median SFH of central galaxies. 
In reality, for a given halo mass the SFR distribution has 
a broad (roughly log- normal) distribution (see Fig. [Ij. 
For a log-normal distribution the average SFR is larger 

than its median value by a factor of e''"^^"''^-' which 
corresponds to an enhancement in the SFDs of ~ O.ldex 
for the typical dispersion, a ~ 0.3, shown in Fig.[TJ Tak- 
ing this into account will increase the predicted SFD, 
especially at z = because of the bi-lognormal distri- 
bution of galaxies in halos of ^ 10^^ /i^^M©. Such an 
increase of ~ O.ldex is insufficient to explain the appar- 
ent discrepancy at high z. An alternative explanation 
might be that we did not include the contribution due to 
satellite galaxies. Albeit small, adding this contribution 
will also slightly increase our SFD model predictions. In 
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Fig. 11. — Instantaneous versus integrated mass ratios. The red and black shaded regions indicate the 68% confidence levels on the 
'integrated' mass ratios between central galaxy and host halo, Mf^c/^hy a-s function of host halo mass, as obtained by Y12 from the SMFl 
and SMF2 data samples, respectively. Different panels correspond to different redshift intervals, as indicated. The solid curve in each panel 
is the 'instantaneous' mass ratio, SFE X /;,, inferred from our model. For comparison, we also show in each panel the model predictions 
assuming "no-evolution", which simply are the Mt^c/M^ (cyan, short-dashed curves) and SFE X (green, long-dashed curves) at ^ = 0.1. 
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Fig. 12. — The cosmic star formation densities. Left-hand panel: comparison of our model predictions with the observational data 
compiled by Ifopkins & Beacom (2006; small dots with error bars). Our model predictions based on SMFl and SMF2 are shown as the 
solid and long-dashed lines, respectively. For comparison, the shaded band indicates the range of SFDs obtained from the 'MIN', 'MAX' 
and 'OBS' models based on SMFl and SMF2, and therefore (roughly) reflects our model uncertainties. Right-hand panel: model predictions 
based on SMFl for the contributions to the total SFD due to halos of different masses, as indicated. 
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addition to these, the high redshift SFD depends signifi- 
cantly on the correction of faint galaxies. As pointed out 
recently by Behroozi et al. (2012) based on some new 
measurements of the cosmic SFD (e.g., Bouwens ct al. 
2012), the data compiled by Hopkins & Bcacom (2006) 
is likely to overestimate the cosmic SFDs in the redshift 
range 3 < z < 8. This discrepancy is largely due to 
the different luminisity cuts in calculating the SFD (see 
Kistler et al. 2009). Taking all these uncertainties/issues 
into account, we conclude that our model seems to pre- 
dict a cosmic SFD somewhat lower than the data at very 
high redshift. The discrepancy can be significantly eased 
either if the stellar mass functions at high redshift have 
a significantly steeper low mass end slope, or if the faint 
end slopes of the galaxy luminosity functions used in the 
SFD measurements arc significantly under-estimated. 

Finally, the right-hand panel of Fig. [12] shows our 
model prediction for how halos of different masses 
contribute to the cosmic SFD. Halos with Mh < 
-|^q10.5 are predicted to only contribute signifi- 

cantly at very high redshifts (z > 10). Halos with masses 
10.5 < \og{Mh/ /i~^M0) < 11.5, on the other hand, are 
predicted to be the main contributors of the cosmic SFD 
at both z ^ and z ~ 7.0. Milky- Way sized halos with 
masses 11.5 < log(M/j/ h~-^MQ) < 12.5 are predicted to 
be the main contributors for most of the history of the 
Universe, all the way from z ~ 7 to z ^ 0.1. Interestingly, 
our model predicts that halos with Mh ^ 10^^'^ /i^^M© 
never contributed significantly (i.e., more than 10%) to 
the cosmic star formation density at any redshift. 

5.6. The fraction of stars formed in situ 

Recent cosmological simulations of galaxy formation 
show that the assembly of massive galaxies, which are al- 
most always ellipticals, consists of "two phases" ; a rapid 
early phase at z > 2, during which stars are formed in 
situ (i.e., within the galaxy) from infalling cold gas, fol- 
lowed by an extended phase during which ex situ stars 
(in the form of satellite galaxies) are accreted (e.g. Oser 
et al. 2010, 2012; Hirschmann et al. 2012). Such a 
two-phase formation scenario for massive galaxies is sup- 
ported by the strong evolution in the observed size-mass 
relation of massive galaxies (e.g., Bezanson et al. 2009, 
and references therein). Our model for the SFHs of cen- 
tral galaxies can predict the fraction of stars formed in 
situ as a function of redshift, halo mass, and stellar mass. 

The total stellar mass of a central galaxy at a given 
redshift zq can be obtained from Eq. (|2ip . The stellar 
mass of stars that formed in situ can be obtained by re- 
placing SFRmax by the real star formation rate SFRqbs , 
so that we can write 

M*,Sf(zo) = / ° SFROBS(^(0) /passivc(tzo ~ t) dt . 
JO 

(27) 

Fig. [13] shows the ratio /msitu = M*,sf(zo)/M*(zo), be- 
tween the mass of stars formed in situ and the total stel- 
lar mass, in both the Alh - z (left-hand panel) and the 
M« - z (right-hand panel) planes. We first focus on the 
most massive galaxies, which are typically ellipticals in 
massive halos. According to our model, at z > 2.5 more 
than 99% of the stars in these galaxies are formed in 
situ. This fraction decreases as a function of redshift, 
dropping to ~ 60% at z = 0. Thus even today's most 



massive ellipticals, for which the accretion of stars from 
(satellite) galaxies is expected to be most important, are 
predicted to have the majority of their stellar mass con- 
tributed by early, in situ star formation (on average). 
For central galaxies in less massive halos, the fraction 
of stars formed in situ is even lower. For a Milky Way 
sized halo of Mh ~ lO^'^ h~^MQ at the present day, more 
than 80% of the stars are expected to have formed in 
situ. Clearly, major mergers of stellar components (e.g. 
major 'dry mergers') of galaxies cannot be a dominant 
mode of stellar mass assembly for galaxies of any stellar 
mass (on average). 

6. SUMMARY 

In this paper we have presented model predictions for 
the star formation histories (SFHs) of central galaxies as 
a function of halo mass. The model is based on self- 
consistent modeling of the conditional stellar mass func- 
tion across cosmic time by Yang et al. (2012). Two key 
ingredients are used in deriving the SFHs. The first is 
the mass assembly histories of central galaxies and their 
accreted satellite galaxies, and second is the local obser- 
vational constraints on the star formation rates of central 
galaxies as function of halo mass. The difference in total 
stellar mass between the accreted and surviving satel- 
lites provides the maximum contribution available to the 
growth of the central galaxy through accretion of stars 
from satellites. The minimum (zero) and maximum con- 
tributions from the accreted satellites correspond to the 
maximum and minimum amounts of stars that formed 
in situ in a central galaxy. As expected, the local obser- 
vational constraints on the star formation rates (SFRs) 
of central galaxy, obtained from the SDSS DR7 group 
catalog, fall between these extrema. 

Using these data, we have obtained median SFHs for 
central galaxies as function of their present-day halo 
mass. We have presented a universal fitting formula that 
adequately describes the dependence of these SFHs on 
halo mass, galaxy stellar mass and redshift. We also used 
this model to make predictions for (i) the star formation 
rates (SFRs), star formation rate densities (SFRDs), gas 
mass depletion rates (GMDRs), and star formation ef- 
ficiencies (SFEs), all as functions of redshift, halo mass 
and stellar mass; (ii) the cosmic star formation rate den- 
sity; and (iii) the fraction of stars that have formed in 
situ over cosmic time (as apposed to have been accreted). 
Our main findings can be summarized as follows. 

1. The SSFR at high z increases rapidly with increas- 
ing redshift [oc (1 -|- z)^'^] for halos of a given mass, 
and slowly with halo mass (ex M°-^^) for a given 
z. This scaling is almost identical to that of the 
specific mass accretion rate of dark matter halos 
(Dekel et al. 2009; McBride et al. 2009; Fakhouri 
et al. 2010), indicating that the SFR of (star- 
forming) central galaxies is largely regulated by the 
rate at which their host halos accrete mass. Such 
a picture has strong theoretical support (e.g., Dut- 
ton, van den Bosch & Dekcl 2010; Bouche et al. 
2010; Dave, Finlator & Oppenhcimer 2012). 

2. The ratio between the SFR in the halo's main pro- 
genitor and the final stellar mass of a galaxy peaks 
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Fig. 13. — Model predictions for the fraction of stars of central galaxies that formed in situ, /insitm as function of halo mass (left-hand 
panel) or stellar mass (right-hand panel) and redshift. Color coding is indicated in vertical bar at the right-hand side. For clarity, the 
purple lines indicate the contours where /insitu = 0.6, 0.7, 0.8, 0.9 and 0.99, respectively. The solid black curves once again show the median 
growth of halo mass (loft-hand panel) and stellar mass (right-hand panel) along the main branch of halos with different present-day masses. 



roughly at a constant value, ^ 10~^-^/i^yr~^, inde- 
pendent of the halo mass and stellar mass of the 
galaxy at the present day. The redshift at which 
this SFR peaks (zpk), however, increases rapidly 
with the present-day halo mass of the galaxy, with 
Zpk ^ 0.5 for Mh = 10" h-^MQ, and Zpk - 3 for 
Mh = IQI^/i-iMq. 

3. More than half of the stars in the present-day 
Universe were formed in halos with masses be- 
tween 10" ! /i-^Mq and lO^^.s /j-iMq in the red- 
shift range 0.4 - 1.9. Halos with masses between 
10"-^ H-^Mq and lO^^-^ H'^Mq dominate the star 
formation rate density of the universe over a large 
range of redshift, from z^lto~5;atz<l the 
star formation rate density is dominated by halos 
with 101°-^ h-^MQ < Mh < 10"-5 H-^Mq; the to- 
tal amounts of stars formed in small halos with 
Mh < 10^°'^ h~^MQ and in massive halos with 
Mh > 10^^-^ H'^Mq are both negligibly small at 
any z < 5. 

4. For individual centrals, the SEE, defined as the 
star formation rate divided by the baryonic ac- 
cretion rate, initially increases, until it reaches a 
maximum, after which it decreases rapidly to a 
quenched state. Eor centrals in present-day ha- 
los with Mh > 10^"^ h~^MQ, quenching occurs 
when their main progenitor halo reaches a mass 
- 10^2 '^ /i-^Mq, at which point the SEE is - 
0.1. Eor centrals in present-day halos with Mh < 
lO^^/i~^M0 the quenching mass shifts to lower halo 
masses at higher peak SEE; at the present, the 
quenching mass is ^ 2 x 10^^ H^^Mq, with a peak 
SEE close to unity. 

5. Whereas the SEE histories of central galaxies 
that end-up in present-day halos with Mh ^ 
10^^ '^ h~^MQ arc consistent with the presence of 



a halo mass floor of ~ 10"'^^ /i~^Mo, as suggested 
by Bouche et al. (2010), our model indicates that 
such a halo mass floor (below which star formation 
is suppressed) needs to be substantially lower, or 
even absent, at high z. 

6. There are some indications that our model may 
underpredict the cosmic star formation density at 
high redshifts (z ^ 3). The discrepancy can be sig- 
nificantly reduced if either the stellar mass func- 
tions at high redshift have a significantly steeper 
low-mass end slope, or the faint-end slopes of the 
luminosity functions used in the SED measure- 
ments arc significantly under-estimated. 

7. At redshift z > 2.5 more than 99% of the stars in 
the progenitors of massive galaxies (mainly ellipti- 
cals) arc formed in situ, and this fraction decreases 
as a function of redshift, dropping to ^ 60% at z = 
0; for a Milky Way sized halo of Mh - lO^^ /j-iMq 
more than 80% of all the stars in the central galaxy 
are formed in situ. Hence, major mergers cannot 
be a dominant mode of stellar mass assembly for 
any stellar mass (on average). 

ACKNOWLEDGEMENTS 

We thank Shiyin Shen, Chunyan Jiang for help in 
dealing with IDL plotting, and Stephane Chariot for 
kindly providing us the passive evolution factors in elec- 
tronic forms. We also thank the anonymous referee 
for helpful comments that improved the presentation of 
this paper. This work is supported by the grants from 
NSEC (Nos. 10925314, 11128306, 11121062, 11233005) 
and CAS/SAEEA International Partnership Program 
for Creative Research Teams (KJCX2-YW-T23). HJM 
would like to acknowledge the support of NSE AST- 
1109354 and NSE AST-0908334. 



20 



Yang ct al. 



O 



8 - 



1 1 1 1 — I — I — r-|; 

M^=10'" h-'Mg. M,= 10"-= h-^Ma 
M^=10"" h-'Mg. M,= 10"° h-=Mg 
M^=10'2« h-'Mj,, M,= 10'°» h-^Mj,— 
_M^=10" = h-'Mg. M,= 10'» h-=Mo 
..M. = 10"'» h-'M„. M.= 10«» h-^M, 




1+Z 




1+Z 



Fig. 14. — Similar to right-hand panel of Fig. |4l except that here the star formation histories are obtained using SMF2, rather than 
SMFl. In the left-hand panel the SFHs are normalized by the stellar masses of the central gal axies at z = 0.1, while the curves in the 
right-hand panel have been normalized by an additional enhancement factor, /cnhi given by Eq. IIA2I I. 
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Fig. 15. — The star formation histories (SFHs) of central galaxies in halos of different present-day masses, as indicated in each panel. 
Here we only show the results obtained from SMF2. In each panel, predictions with MIN and MAX assumptions for the star formation 
rate are bridged with shaded areas (see text for more details). In panels for halos with Af;, > lO"*^^ h~^MiT), local observational constraints 
are shown as the vertical shaded histograms (distributions) and stars (median). The dotted, dashed, and long-dashed lines are the MIN, 
MAX and OBS fits to the SFHs discussed in the text. 



APPENDIX 

A. THE SFH FITTING FORMULA FOR SMF2 

The discussion in the main text is mainly based on SMFl. As we have seen, the SFHs obtained from SMFl and 
SMF2 have systematic differences (see, e.g., Fig. [2]). For completeness, this Appendix presents our model for the SFHs 
of central galaxies based on SMF2, rather than SMFl. Let us first look at the amplitude of the SFHs, shown in the 
left panel of Fig. [TH Compared to the results obtained from SMFl (sec the right panel of Fig. 2]), the peak values 
of the SFRs obtained from SMF2 have slightly larger variation between different halos masses. Especially the very 
massive galaxies in SMF2 have significantly higher SFH peaks than those in SMFl. In order to properly model the 
SFH peaks in the SMF2, we introduce a stellar mass dependent enhancement factor /cnh, so that we can properly 
model the peak amplitudes. By fitting to the SFH peaks obtained from SMF2, we get the following relation between 
SFRpk and M^^o.i- 

SFRpk _ M*.o.i tnj \ 

[Moyr-i] " 109-4 /i-2Mq^°'^'^^^^*'°-^^ ' ^^^^ 
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where 

/o„h(M*,o.i) = 1 + (A/,,o.i/10"-2 h-^M^f , (A2) 

is a stehar mass dependent enhancement factor. The performance of the fitting results is shown in the right panel of 
Fig. [TH As one can see, the model describes the amplitudes remarkably well. 
The shape of the SFH obtained from SMF2 can still be modeled using the same form as Eq. (|18|) , 

SFR(A4, z) = SFRp, X [~^^^^^^^/^^^] , (A3) 

and here with (T(zpk) given as follows. For z > Zpk, 

(7(zpk) = 0.146(1 + Zpk)-°-"^ (A4) 

while for z < Zpk it is 

f 0.0857(1 + zpk)°-^^^ (OBS) 
a(zpk) = <^ 0.168(1 + Zpk)-°-208 (MAX). (A5) 

[ 0.346(1 + zpk)^^-^^ (MIN) 

The predictions of this fitting model are shown as the long-dashed (OBS), dashed (MAX) and dotted (MIN) curves 
in Fig. [13] in comparison with the results obtained directly from SMF2. 
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